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Abstract 

We study exact renormalisation group equations for the 3d Ising universality class. 
At the Wilson-Fisher fixed point, symmetric and antisymmetric correction-to-scaling 
exponents are computed with high accuracy for an optimised cutoff to leading order 
in the derivative expansion. Further results are derived for other cutoffs including 
smooth, sharp and background field cutoffs. An estimate for higher order corrections 
is given as well. We establish that the leading antisymmetric corrections to scaling 
^ \ are strongly subleading compared to the leading symmetric ones. 
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1. Introduction 



Many physical systems with short range interactions and a scalar order parameter 
display Ising universal behaviour close to the critical point. Initially introduced for the 
study of magnetic systems, the Ising model also describes the physics of the liquid-gas 
phase transition, transitions in binary mixtures and in Coulombic systems In high 
energy physics, Ising universal behaviour is expected in various theories including the QCD 
phase transition with finite quark masses 2], the chiral phase transition of QCD 0,0, and 
the electro-weak phase transition 



The original Ising model has a global Z 2 symmetry. However, many systems in the Ising 
universality class like the liquid-gas and the electro-weak phase transition do not possess 
the Z 2 symmetry away from the critical point. The presence of operators odd under Z 2 lead 
to additional corrections-to-scaling exponents. In principle, deviations from Z 2 symmetric 
scaling are detectable experimentally. Antisymmetric corrections to scaling ~ L -0,5 have 
been revealed in a Monte-Carlo simulation of the electro- weak phase transition . Previous 
theoretical studies of antisymmetric corrections to scaling are based on the e-expansion 0, 
the scaling field method p, and the Wegner-Houghton equation p|. 

In this Letter, we study corrections to scaling for the 3d Ising universality class using 
the Exact Renormalisation Group, which is based on the Wilsonian idea of successively 
integrating out momentum modes (see 0] for reviews and references therein). This 
approach is implemented through an infra-red cutoff which, within a few constraints, 
can be chosen freely. The strengths of the method are its flexibility and its numerical 
stability. Furthermore, a general optimisation procedure is available, enabling a choice 
of the infra-red cutoff best suited for the problem at hand 0. To leading order in a 
derivative expansion, we employ an optimised cutoff and compute the first six subleading 
corrections-to-scaling exponents with high accuracy. We also obtain results for smooth, 
sharp and background field cutoffs, and estimate higher order corrections. 



2. Renormalisation group and critical exponents 

Renormalisation group methods have been used very successfully in the computation of 
universal observables at second order phase transitions. A particularly useful approach is 
the Exact Renormalisation Group (ERG), based on the Wilsonian idea of integrating out 
momentum modes within a path integral representation of quantum field theory y|. In its 
modern form, the ERG flow for an effective action for bosonic fields if is given by the 
simple one-loop expression 

d t T k [ip] = ±Tr(T^ + Ry 1 d t R. (1) 

Here, t = In k is the logarithmic scale parameter, the trace denotes a momentum trace and a 
sum over indices, r( 2 )[<£>](p, q) = 5 2 T /5cp(p)5cp(q), and R is an infra-red momentum cutoff at 
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the momentum scale k. The flow (JIJ interpolates between an initial (microscopic) action in 
the ultra-violet and the full quantum effective action in the infra-red. At every momentum 
scale k, ((J) receives its main contributions for momenta about p 2 ~ k 2 . The regulator can 
be chosen freely and allows for an optimisation of the flow within general approximations 
jj. The optimisation entails improved convergence and stability properties of the flow. 
In combination with the numerical stability of the flow, it provides the basis for reliable 
predictions based on systematic approximations to ((TJ). An important non-perturbative 
approximation scheme is the derivative expansion [10| . To leading order, the local potential 
approximation consists in the Ansatz 

T k = J d d x (u k (ip) + ^dpipd^ (2) 

for the effective action. It implies that higher order corrections proportional to the anoma- 
lous dimension 77 of the fields are neglected. For the Ising universality class, rj is of the order 
of a few percent. Inserting the Ansatz (j2J) into the flow equation (0) and evaluating it for 
constant fields leads to the flow for the effective potential Uk- We introduce dimensionless 
variables u(4>) = U k /k d and <p — ^pk 1 ~ d ^ 2 . Then, finding a fixed point amounts to solving 
dtu = 0. To that end we employ a polynomial truncation of the fixed point potential, 
retaining vertex functions n up to a maximum number n trunc , 

^trunc 

= E (3) 

71=1 

The potential has been normalised as u(<f) = 0) = 0. The Ansatz (jSJ) leads to n trunc ordinary 
differential equations dtTi = /3j. In three Euclidean dimensions, the flow equation exhibits the 
nontrivial Wilson-Fisher fixed point 7^ 0. Universal critical exponents and corrections- 
to-scaling exponents are obtained as the eigenvalues of the stability matrix at criticality 
Mij = d{3i/dTj\*. For convenience, we introduce the set of 0-even couplings A n = r 2n and 0- 
odd couplings A n = r 2n _i. Under reflection in field space the couplings and their /9-functions 
behave as 

The scaling solution is symmetric under — > —0. Hence all 0-odd couplings A* vanish 
at the fixed point. The computation of critical exponents is simplified by observing that 
/%(A,A) = — P\(X, — A) for all A. This follows from (JHJ) and (jlj). In particular, f3\ vanishes 
identically at A = 0, where all derivatives of (3\ w.r.t. the symmetric couplings vanish at 
the fixed point. Therefore, the stability matrix M at criticality simplifies and becomes 
equivalent to the matrix 

(AB\ 

loci 



(5) 



with A = dfix/dXl*, B = dfix/dXl* and C = dflx/dX^. In consequence, the eigenvalues 
of M reduce to those of the sub-matrices A and C. The matrix A carries the information 
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Figure 1: The exponents u, u and 002 from a polynomial 
truncation to order n trunc about vanishing field. 



about the critical exponent v and the symmetric corrections to scaling, while the matrix C 
contains the information about antisymmetric corrections to scaling. 



3. Results 



In this section, we present our results for the universal eigenvalues of the stability matrix. 
The fixed point is determined in truncations up to n trunc = 40. The optimised regulator 
Rp ut = (k 2 — q 2 )6{k 2 — q 2 ) is employed to improve the convergence and stability of the flow 
The stability matrix is evaluated with two different methods: an expansion in 



powers of the field about vanishing field, and an expansion in Legendre polynomials. The 
latter case involves an integration in field space. 

Our numerical results are given in Tab. 1 and 2, and in Figs. 1-4. First, we discuss our 
results for the 0-even corrections to scaling. In Fig. 1, the three leading eigenvalues of A are 
displayed, v" 1 , 00 and U2. Notice that the eigenvalues are identical for truncations (n, n + 1), 
if n is even. The reason for this is simple: increasing the truncation by a 0-odd coupling does 
neither change the dimension of the matrix A, nor the numerical values of the fixed couplings 
(because 0-odd couplings vanish at the fixed point). Hence, the 0-even eigenvalues remain 
unchanged, as is clearly seen in Fig. 1. The numerical values for the 0-even eigenvalues are 
identical to those which are found in a polynomial expansion in p = <j) 2 /2. 
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Figure 2: The exponent oj\ (see Fig.l). 



0-even 


0-odd 






Vh 


-2.5 


V 


0.6495 


2/shift 


-0.5 


UJl 


0.655 


Ui 


1.88 


^2 


3.18 




4.5 


U 3 


5.9 


u 3 


7 



Table 1: 0-even and 0-odd eigenvalues (-R op t, ^trunc — 

40). 

Next, we turn to the 0-odd corrections to scaling (Figs. 2-4). We have computed the 
eigenvalues of the matrix C for n trunc up to n trunc = 40. We find two eigenvalues = —5/2 
and 2/ s hift = —1/2, related to redundant operators i6|. All other eigenvalues are positive. We 
denote them as u n . The leading non-trivial 0-odd correction-to-scaling exponent Uj\ as a 
function of the truncation is displayed in Fig. 2. (In the literature, U)\ is sometimes denoted 
as uja or u 5 .) Our results for oj 2 and lj 3 are given in Fig. 3 and 4, respectively. Notice that 
the pattern of the results, with increasing truncation, is similar to the pattern found in the 0- 
even sector. The results for two subsequent truncations (n, n + 1) for n odd are close to each 
other for sufficiently large truncation. The reason for this is the following: increasing the 
truncation by a 0-even coupling does not change the dimension of the matrix C. However, 
it does change the numerical value of the fixed point, and hence the eigenvalues of C. With 
increasing truncation, the numerical change within the 0-even couplings at the fixed point 
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Figure 3: The exponent ui2 (see Fig.l). 
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Figure 4: The exponent 0)3 (see Fig.l). 



is very small and eventually, the eigenvalues of C become insensitive to the addition of a 
0-even coupling. This is nicely observed in the results presented in Figs. 2-4. Comparing 
the symmetric with the antisymmetric corrections to scaling, the general pattern we find is 
that < iO\ < uji < ■ ■ ■ < u n < uj n < ■ ■ -. 



0-even 


0-odd 






Vh 


- 2.5 


V 


0.649562 


?/shift 


-0.5 


lO\ 


0.6557 


Qi 


1.886 


U2 


3.180 


C02 


4.524 


UJ 3 


5.912 


UJ 3 


7.33 



Table 2: c/>-even and </>-odd eigenvalues (i? pt) ^trunc — 

22, 



0.46). 



We have also computed the critical exponents by using the approach of [7(. Expanding 
the scaling potential and the eigenperturbations in terms of orthogonal polynomials 
(Legendre polynomials) implies that the matrix elements of © involve an integration in 
field space 4> G [— <fi max _, <fi m ax\- The weak dependence on max is fixed by requirin g th at the 
0-even eigenvalues agree to high accuracy with the known values obtained in |l2| using 
an expansion in p = <j) 2 /2 about the potential minimum p = p . This procedure improves 
the numerical convergence. Our results for the eigenvalues are given in Tab. 2. They are 
consistent with and have a higher accuracy than those given in Tab. 1. 
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Figure 5: The exponent Co\ from a polynomial expansion up to 
order ntruno and in comparison for the sharp, the background 
field (bg), the optimised (opt) and the quartic cutoff (see text). 



Next we discuss our results based on other regularisations including the power law cutoff, 
the sharp cutoff and a background field cutoff. Varying the momentum dependence of the 
regulator R from "smooth" to "sharp" allows for an estimate of higher order corrections 
due to operators neglected in the present approximation, e.g. (T^f. We have computed the 
0-even and 0-odd eigenvalues for a smooth power-like regulator -R qU artic = k 2 ■ (q 4 /k 4 ), for 
the sharp cutoff -R s harp — nm a^oo a k 2 Q(k 2 — q 2 ), and for a background field cutoff i?b g - 
Results are summarised in Fig. 5 and Tab. 3. 



The power-law cutoff -R qU artic is optimised in the sense of |9|. We find that the flow 
based on -R qua rtic has similar stability and convergence properties as the flow with -R op t, 
e.g. Fig. 5. Also, the numerical results as given in Tab. 3 are within 5% or less to each 
other. The sharp cutoff does not lead to an optimised flow It displays instabilities 
within a local polynomial expansion about vanishing field 3, 1A \ ■ This is confirmed in our 



analysis. In Fig. 5, the 0-odd eigenvalue uj\ is displayed up to a truncation n trunc = 30. The 

local field expansion based on -R s har P oscillates in the eight-fold pattern (+ + + H ) 

about its mean value, reminiscent of the four-fold pattern in an expansion in 2 , 0]. The 
expansion fails to converge at the present order. The asymptotic value for uj\ is indicated 
by the dashed line. Moreover, the eigenvalues for ri trU nc = 16, 17, 24 and 25 possess a small 
imaginary part, which is not displayed in Fig. 5. These findings are a reflection of an 
intrinsic instability of the sharp cutoff flow. Furthermore, the sharp cutoff value for the 



jp c 

leading critical exponent v is ~ 8% larger than the value for R opt jllj . and ~ 10% larger 
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than the physical value. These properties indicate that quantitative results from the sharp 
cutoff flow within a given truncation, although qualitatively correct, are less reliable than 
those by optimised cutoffs. 

The instability in the stability matrix of the sharp cutoff flow is removed by expanding 
the fixed point potential and the eigenperturbations in terms of Legendre polynomials. For 
the eigenvalues, we adjust max as described above to improve the numerical convergence. 
Our results are given in Tab. 3, and by the dashed line in Fig. 5. In the 0-even sector, our 
results agree to all significant digits with those by Morris (quartic cutoff) [15| . Cornelias and 
Travesset (sharp cutoff) 2| and Litim (optimised cutoff) This provides a non-trivial 
consistency check, because the numerical methods employed in 15, 3|, in [l^l, and here, 
are all different. The main new results concern the eigenvalues in the 0-odd sector, where 
we also confirm the first two eigenvalues by Tsypin (sharp cutoff) 0| (see also [l3|)- 



0-even 




-Rquartic 


-Ropt 




0-odd 


-Rsharp 


-Rquartic 


-Ropt 


-R bg 


V 


0.6895 


0.6604 


0.649562 


0.625979 












Ui 


0.5952 


0.6285 


0.6557 


0.762204 


Qi 


1.691 


1.812 


1.886 


2.163 


UJ2 


2.838 


3.048 


3.180 


3.6845 




3.998 


4.32 


4.524 


5.313 


LU 3 


5.18 


5.63 


5.912 


7.038 


UJ 3 


6.38 


6.96 


7.33 


8.85 



Table 3: 0-even and 0-odd eigenvalues for the sharp, quartic, optimised and background 
field cutoff, using max = 0.43,0.45, 0.46, and 0.46, respectively (see text). 



Now we proceed with the background field flow, where plain momenta q 2 in the regulator 
are replaced by_T( 2 )[0], the full inverse propagator evaluated at some background field 
jl8^ (see also jl^). Identifying the background field with the physical mean leads to a 
partial diagonalisation, which should further stabilise the flow. Background field flows 
are closely related to the proper-time renormalisation group of Liao [20], to which they 
reduce once an additional flow term proportional to c^r^ 2 -* is dropped. Implicit to this 
approach is that differences between fluctuation and background field are neglected — an 
approximation, which in the present theory becomes exact in the infra-red limit. Hence, 
as detailed in 18(, we expect this approximation to be viable in the vicinity of a critical point. 

Here, we use the flow d t Tk = Tr exp —T^/k 2 to leading order in a derivative expansion. 
It is a background field flow in the proper-time approximation with regulator i?bg given 
by (13), (14) of Amongst the proper-time flows, it has best stability properties [2]J 
(see also [Hj). This is reflected by the very fast convergence of ui\ with the truncation 
(Fig. 5). We also stress that the first two eigenvalues in the 0-even sector, which agree 
with earlier results in j^, 2^|, are very close to the physical values. The further subleading 
corrections-to-scaling exponents are increasingly larger than the values for R ov t- This trend 
is indicative for the potential effect of higher order corrections. 
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Finally we comment on the Polchinski renormalisation group |23j]. It is related to the 
flow (JIJ by a Legendre transform and additional field rescalings. In consequence, both 
methods have inequivalent derivative exp ansions. To leading order, the Polchinski flow 
is independent of the regularisation [24j |. in contradistinction to the present approach, 
e.g. [12] . For R opt , critical exponents in the 0-even sector agree to high precision with those 
from the Polchinski flow. The numerical value for the 0-odd eigenvalue a)i given in Tab. 2 
for -R op t) also agrees with preliminary results from the Polchinski flow [25j. If these find- 
ings persist, they confirm the deeper link between the two methods even for the 0-odd sector. 



4. Discussion and conclusion 

We have studied symmetric and antisymmetric corrections to scaling at criticality for 
systems belonging to the 3d Ising universality class. The first six subleading universal 
corrections-to-scaling exponents are obtained from an exact renormalisation group. Best 
results are achieved for optimised flows, which have enhanced convergence and stability 
properties. In addition, we have assessed the cutoff dependence for smooth, sharp and a 
background field cutoff. This study also served as an indicator for higher order effects. 
Results from the standard and the background field flows have to be seen on slightly 
different footings due to qualitative differences in the approximations. 

For the optimised flow, the leading symmetric and antisymmetric correction-to-scaling 
exponents are uj = 0.6557 and uja = 1.886. For different regularisations ranging from sharp 
to optimised cutoffs and including (excluding) the background field flow, the exponents 
vary between uj « 0.60 - 0.76 (0.60 - 0.66) and uj a « 1.7 - 2.2 (1.7 - 1.9). Higher order 
corrections due to a non-vanishing anomalous dimension lead to w ~ 0.8, and similar 
corrections are expected for uja- Expressed in terms of the exponent = uja v, our results 
1.22 for the optimised, A^ « 1.2 for the quartic, A^ « 1.17 for the sharp and 
Aa ~ 1-35 for the background field cutoff. This compares well with A^ rs 1.3 which is often 
assumed in the analysis of experimental data, e.g. |26J]. The leading symmetric corrections 
to scaling are A = ujv « 0.42 — 0.48, increasing towards A « 0.52 once anomalous 
dimensions are taken into account. This cutoff dependence indicates the expected size of 
higher order effects. Our results for uja are consistent with the estimate uja > 1.5 based on 
3 ade resummation of the e-expansion |5j, and with uja ~ 2.4 from the scaling field method 
(J. We notice that all sharp cutoff eigenvalues are systematically smaller than those from 
any other cutoff. This reflects, we believe, the notoriously poor convergence behaviour of 
sharp cutoff flows. 

In conclusion, we have established that the leading antisymmetric corrections to scaling 
are consistently suppressed compared to the leading symmetric ones. Within the errors, 
the exponent uja is more than twice as big as uj. Hence, the scaling behaviour ~ L~ a5 as 
seen in a Monte-Carlo simulation of the electro-weak theory clearly dominates over both the 
leading symmetric ~ L~ u and antisymmetric ~ L^ U1A corrections to scaling and therefore 
cannot be explained with the exponent uja- 
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